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We show that spatially local, yet low-energy, fluctuations can play an essential role in the physics 
of strongly correlated electron systems tuned to a quantum critical point. A detailed microscopic 
analysis of the Kondo lattice model is carried out within an extended dynamical mean-field approach. 
The correlation functions for the lattice model are calculated through a self-consistent Bose-Fermi 
Kondo problem, in which a local moment is coupled both to a fermionic bath and to a bosonic 
bath (a fluctuating magnetic field). A renormalization-group treatment of this impurity problem — 
perturbative in e = 1 — 7, where 7 is an exponent characterizing the spectrum of the bosonic 
bath — shows that competition between the two couplings can drive the local-moment fluctuations 
critical. As a result, two distinct types of quantum critical point emerge in the Kondo lattice, 
one being of the usual spin-density-wave type, the other "locally critical." Near the locally critical 
point, the dynamical spin susceptibility exhibits lu/T scaling with a fractional exponent. While the 
spin-density-wave critical point is Gaussian, the locally critical point is an interacting fixed point at 
which long-wavelength and spatially local critical modes coexist. A Ginzburg-Landau description for 
the locally critical point is discussed. It is argued that these results are robust, that local criticality 
provides a natural description of the quantum critical behavior seen in a number of heavy-fermion 
metals, and that this picture may also be relevant to other strongly correlated metals. 
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I. INTRODUCTION 

Non-Fermi-liquid properties have been seen experi- 
mentally in a number of strongly correlated electron 
systems. 1-3 They pose fundamental questions about how 
electron correlations lead to new electronic states of mat- 
ter and new elementary excitations. One mechanism for 
the breakdown of Fermi-liquid theory is quantum criti- 
cality. While the extent to which an underlying quantum 
critical point (QCP) plays a role in the non-Fermi-liquid 
behavior of high-temperature superconductors remains 
a subject of debate, the situation is much clearer in 
heavy-fermion metals. Among the many heavy-fermion 
materials in which non-Fermi-liquid properties have been 
seen, 1,5-15 magnetic QCPs have been explicitly identi- 
fied in a number of stoichiometric or nearly stoichio- 
metric materials 5-8 including CeCue-z Au^, CePd2Si2, 
Celri3, and YbRli2Si2. In each of these systems, the zero- 
temperature magnetic ordering transition appears to be 
continuous. Non-Fermi-liquid behavior, usually seen in 
the temperature dependences of transport and thermo- 
dynamic properties, arises near the QCP. For instance, 
the resistivity is linear (or close to being linear) in T, 
and the temperature dependence of the specific heat co- 
efficient C /T is either logarithmically singular or nonan- 
alytic with a finite zero-temperature limit. Away from 
the quantum critical regime, there is a gradual recovery 
of Fermi-liquid behavior (albeit still with a large effec- 
tive mass). These materials provide controlled settings 
in which to study not only basic issues concerning mag- 
netism in heavy fermions, but also the physics of QCPs 
in strongly correlated metals in general. 

Inelastic neutron scattering directly probes the critical 
fluctuations at heavy-fermion magnetic QCPs. Partic- 



ularly detailed measurements have been performed on 
CeCug-ajAua;, which can be tuned from a paramagnetic 
metal to an antiferromagnetic metal by varying x and/or 
an applied pressure. The results 16 " 19 are striking in a 
number of ways: (1) The frequency dependence of the 
dynamical spin susceptibility displays a fractional expo- 
nent. (2) The same exponent also appears in the tem- 
perature dependence; the dynamical spin susceptibility 
has ui/T scaling. (3) This fractional exponent describes 
the variation with frequency and temperature not only at 
the ordering wavevector, but essentially everywhere else 
in the Brillouin zone. 

The existence of the fractional exponent is a surprise. 
It is standard to assume 20-25 that the critical theory for 
magnetic QCPs in metals is a cf> 4 theory, describing the 
long-wavelength fluctuations of the magnetic order pa- 
rameter. Landau damping makes the dynamic exponent 
z larger than 1. In the antiferromagnetic case, z = 2, so 
the effective dimensionality of the <j> theory, d e g = d + z, 
is greater than (equal to) its upper critical dimension for 
spatial dimensionality d — 3 (d — 2). Consequently, in 
the standard picture all the nonlinear couplings are irrel- 
evant in the renormalization-group (RG) sense, the fixed 
point is Gaussian, and the frequency exponent must take 
its mean- field value of 1. 

Likewise, the existence of u/T scaling is surprising. 
For a Gaussian fixed point, the spin relaxation rate is 
determined by nonlinear couplings that are irrelevant in 
the RG sense. This implies that the relaxation rate is 
superlinear in temperature, 26 whereas uj/T scaling can 
arise only for a linear relaxation rate. 

Both of these features, then, imply that the critical 
point in CeCue-^Aua; has to be an interacting one. What 
is the new physics that is responsible for the interact- 
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ing part of the critical theory? A clue is provided by 
the third feature of the experimental data mentioned 
above. The fact that the fractional exponent in the fre- 
quency/temperature dependence occurs also at generic 
wavevectors (far away from the ordering wavevector) sug- 
gests that the origin of the fractional exponent lies in 
some type of local physics. 17 

Hcavy-fermion metals undergoing a magnetic quan- 
tum phase transition should be well-described by the 
Kondo lattice model of local moments interacting with 
conduction-band electrons. It is natural, therefore, to 
suspect that the local physics responsible for the anoma- 
lous properties of CeCu6-xAu x involves the magnetic 
moments. In Ref. 27, we briefly reported an interacting 
critical point in Kondo lattice systems. (See also Ref. 28 
for more qualitative discussions.) Our analysis was car- 
ried out using an extended dynamical mean field theory 
(EDMFT), supplemented by Ginzburg-Landau consider- 
ations. Other approaches, based on a large-iV formula- 
tion, have also been developed to study this problem, 29 ' 30 
although they have yet to yield any new type of QCPs. 
A scenario of three-dimensional electrons coupled to two- 
dimensional Gaussian spin fluctuations 19 ' 31 ' 6 has also 
been proposed, but it does not address the local and other 
non-Gaussian aspects of the critical dynamics. 

The present paper describes in detail our microscopic 
analysis of the Kondo lattice. The interplay of local and 
spatially extended physics is treated using the EDMFT 
developed in Refs. 32-34. The central result is the find- 
ing of two types of quantum phase transition. The first, 
corresponding to the Gaussian picture, describes a tran- 
sition of the usual spin-density-wave (SDW) type, and is 
found to occur when the underlying spin fluctuations are 
three-dimensional in character. We dub the other type 
of transition "locally critical," because long-wavelength 
and spatially local critical modes coexist at the QCP. 
This type of transition occurs if the spin fluctuations are 
(quasi-) two-dimensional. Near the locally critical point, 
the dynamical spin susceptibility exhibits lu/T scaling 
with a fractional exponent. By analyzing the general 
form of the interacting critical theory, we argue that the 
existence of two types of critical point is valid beyond 
the approximations contained in our microscopic calcu- 
lations. The locally critical picture provides a natural 
description of the quantum critical behavior seen in a 
number of heavy-fermion metals, and is argued to be of 
importance in the broader context of strongly correlated 
electron physics. 

The remainder of the paper is organized as follows. 
We introduce the model and its EDMFT formulation 
in Section II. The EDMFT equations have the con- 
tent of a Bose-Fermi Kondo model supplemented by self- 
consistency conditions. In Section III, we present a de- 
tailed RG analysis of the Bose-Fermi Kondo model, when 
the bath spectral functions are assumed to take power- 
law forms. A critical point is identified in the impurity 
model, and the spin correlation functions near this criti- 
cal point are calculated. Using this solution to the impu- 



rity model, we study in Section IV the self-consistent 
EDMFT problem. Two types of QCP are identified, 
and their zero-temperature dynamics are determined. 
The finite-temperature spin dynamics are studied in Sec- 
tion V. A Ginzburg-Landau analysis is presented in Sec- 
tion VI. Section VII compares our theoretical results 
with experiments in heavy-fermion systems. Section VIII 
contains concluding remarks and discusses the general 
relevance of critical local physics. Details of the deriva- 
tion of the RG equations, of the determination of the sep- 
aratrix between the two stable phases of the Bose-Fermi 
Kondo model, and of the calculation of correlation func- 
tions at the critical point of the impurity model are given 
in Appendices A, B, and C, respectively. Finally, Ap- 
pendix D presents the derivation of correlation functions 
for the long-ranged one-dimensional spherical model with 
an interaction that decays with distance either as a pure 
power-law or with logarithmic corrections. 

II. THE MODEL AND FORMALISM 

A. Kondo Lattice Model 

The Kondo lattice model is specified by the Hamilto- 
nian 

ija i ij 

At each lattice site i is located a spin-i local moment 
S i; which interacts on-site with the spin s C:i of the 
conduction (c) electrons. The tight-binding parameters 
Uj determine the dispersion ej. and, hence, the bare 
conduction-band density of states 

Po (e) = J2^-e k ). (2) 

k 

We will treat only cases where the average number of 
conduction electrons is less than one per site. All the 
phases described below are metallic. 

Two processes oppose each other in the Kondo lat- 
tice. First, each magnetic moment couples locally with 
strength Jk to the spins of the conduction electrons. 
We consider only positive (antiferromagnetic) Jk val- 
ues, which can lead to Kondo-quenching of the local 
moments. Second, the magnetic moments interact with 
each other through the nonlocal, RKKY coupling 7^ . In 
Eq. (1), we have introduced an explicit RKKY interac- 
tion term in addition to that induced implicitly by the 
Kondo coupling Jk- This is done so that the dynamics 
associated with the RKKY interactions can be incorpo- 
rated into the EDMFT. 32 " 34 (The EDMFT ensures that 
there is no double-counting of the implicit and explicit 
RKKY interactions. 32 ) For the systems under considera- 
tion, / q , the spatial Fourier transform of Jy, is most nega- 
tive when q is equal to the antiferromagnetic wavevector, 
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The competition between the Kondo effect and the 
RKKY interaction is expected to iead to a magnetic 
quantum phase transition. 35,36 Traditionally, it is as- 
sumed that the local moments are quenched, not only on 
the paramagnetic side but also through the transition. In 
this picture, the quantum phase transition corresponds 
to an SDW instability of heavy quasiparticles produced 
by the quenching, and the quantum critical behavior is 
necessarily that of a Gaussian fixed point. By contrast, 
our analysis makes no prior assumption that Kondo res- 
onances are fully developed. It should be stressed — as a 
discussion at the end of Section IV A will make clear — 
that the critical point is always reached by tuning just 
one parameter: the ratio of the effective RKKY interac- 
tion and an effective Kondo coupling. 



B. Extended dynamical mean field theory 

In the EDMFT, 32 ~ 34 all the correlation functions of 
the Kondo lattice can be calculated in terms of a self- 
consistent Kondo impurity model. The EDMFT general- 
izes the standard dynamical mean-field theory 37 in that it 
incorporates the quantum fluctuations produced by any 
intersite interaction; in the Kondo lattice model (1) this 
is the RKKY interaction. A key quantity introduced in 
the EDMFT is the "spin self-energy" M(uS), in terms of 
which the momentum-dependent dynamical spin suscep- 
tibility can be written 



x(q,w) 



1 



M{u) 



(3) 



(More specifically, the spin self-energy is defined in terms 
of an effective spin cumulant that is /-irreducible. 32 ) The 
conduction-electron Green's function is, as usual, ex- 
pressed in terms of the conduction-electron self-energy 
E( W ): 



G(k,w) 



H - e k - 



(4) 



A crucial advantage of the EDMFT is that, unlike stan- 
dard RPA approximations, the damping of spin fluc- 
tuations is not restricted to take place via decay into 
quasiparticle-quasihole pairs. Instead the (dynamical) 
solution of the problem will dictate the nature of the low- 
energy many-body excitations and the associated form 
of the spin damping. At the same time, it is impor- 
tant to bear in mind the central approximation made in 
the EDMFT, namely, both £ and M are taken to be 
momentum- independent . 

Details of the formalism have been published in 
Ref. 32. Here we stress that the EDMFT can be adopted 
as a conserving resummation of diagrams for finite- 
dimensional systems. As discussed in Ref. 32, the dia- 
grams retained in the EDMFT form an infinite series, 
which is illustrated in Fig. 1. (Note that standard Feyn- 
man diagrams can be used for the Kondo lattice model 
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FIG. 1. Single-site, two-site, and three-site diagrams for 
the Luttinger-Ward potential in the extended dynamical 
mean-field theory. A dashed line represents an intersite in- 
teraction; this is Iij for the Hamiltonian given in Eq. (1). A 
blob contains all the on-site diagrams, with the associated 
(fully dressed) fermion Green's function. The series extends 
to diagrams involving an infinite number of sites. 



once the spin of the local moment is represented in terms 
of pseudo-fermions without a constraint; one such repre- 
sentation suitable for the spin-i case we are considering 
is the Popov representation. 38 ) 

The spatial dimensionality d and other aspects of the 
lattice structure are encoded in the form of the tight- 
binding parameter t ij and the RKKY interaction 7y . In 
particular, a quantity that will play a crucial role in our 
analysis is the "RKKY density of states," 



(5) 



For any finite-dimensional system, I q is nonzero only 
over a finite range. In other words, the support of pi(e) 
is bounded. This ensures that there is a finite region of 
parameter space over which x(q, to) is positive and non- 
singular; hence, a stable paramagnetic solution exists. 
An alternative construction of the EDMFT, also given 
in Ref. 32, scales the intersite parameters such that the 
nearest- neighbor coupling goes as Iuj\ cx 1/ \fd, and then 
takes the limit d ~ * oo. In this case, a special choice of 
lattice has to be made in order to obtain a stable para- 
magnetic solution. 



C. Extended dynamical mean- field equations for the 
Kondo lattice 

The EDMFT equations for the Kondo lattice model 
can be expressed in terms of an effective action for a 
single lattice site (the "impurity" site): 



^imp — ^top 







dr Jk S • s c 
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rP fP 
dr dr' 
Jo Jo 



+ ^(t)-Xo 1 (t-tI)S(t>) 



(6) 



where (3 = l/fc^T, <S t o P describes the Berry phase of the 
local moment, and Gq 1 and x^ 1 are Weiss fields. 

Equivalently, the effective impurity action can be rep- 
resented in terms of the following effective impurity 
Hamiltonian: 

Wimp = Jk S • S c + E p C pa C pa 
p.a 

w p 4>£ ■ <pp. (7) 

p p 

This is the Bose-Fcrmi Kondo Hamiltonian, describing a 
local moment coupled not only to the spin of a fermionic 
bath (c pCT ), as in the usual Kondo problem, but also to 
a dissipative vector-bosonic bath (<j) p ). Physically, the 
vector bosons describe a fluctuating magnetic field gen- 
erated via the RKKY interaction by the local moments 
at all other lattice sites. The dispersions, E p and w p , as 
well as the coupling constant g, are such that integrat- 
ing out the two baths in the Hamiltonian representation 
yields the effective action (6). Specifically, 



* 2 E 



(iVm) 2 - Wl 



i 



(8) 



= G {iuj n ), 



3. Determine the parameters in the trial forms 
of Xq {i y m) and Gq (iu) n ) by demanding self- 
consistency. The self-consistency conditions 



XimpH 



G{ mp (ui) 



q 



k 



x(q,w), 



G(k,w), 



(13) 



amount to the requirement that each local correla- 
tion function is equal to the wave-vector average 
of the corresponding lattice correlation function. 
The lattice spin susceptibility and lattice Green's 
function, xfawm) an d G(k, iu> n ), respectively, are 
given in Eqs. (3) and (4). 

4. Once self-consistency is achieved, identify the im- 
purity quantities Ximp and G; mp with the corre- 
sponding local (on-site) quantities xioc and G\ oc of 
the lattice problem. Likewise, identify M; mp and 
Si m p with the corresponding lattice quantities M 
and S, and hence use Eqs. (3) and (4) to calculate 
the lattice dynamical spin susceptibility x(q, iv m ) 
and the lattice Green function G(k, iw n ). 

The derivation of Eqs. (3), (4), and (11) has been given 
in Rcf. 32. For notational simplicity, in the remainder 
of the paper we will not differentiate between impurity 
quantities and local quantities. In particular, the impu- 
rity spin-spin correlation function determined from the 
Bose- Fermi Kondo model will simply be called xioc- 



where v m and u> n are bosonic and fermionic Matsubara 
frequencies, respectively. In the following, we will inter- 
change between the two representations based on nota- 
tional convenience. 

The self-consistent procedure goes as follows: 

1. Put in trial forms for Xq 1 (w m ) and Gq 1 (uj„). 

2. Solve the impurity Kondo problem to determine the 
impurity correlation functions, 



and 



Ximp(r) = (T t S*(t)S*(0)> 



Gimp(r) = -(T tCct (t)4(0)). 



(9) 



(10) 



This step also determines the spin self-energy 
through a Dyson-like equation for the impurity 
problem, 

M imp (iv m ) = Xo l {ivm) H j-. y (11) 

Ximp \}^vn ) 

and a conduction-electron self-energy 

Simp(^n) = G^ilUn) - — \-. -. (12) 

Lrimpif^n ) 



III. SOLUTION OF THE IMPURITY PROBLEM 

Our goal is to determine the universal low-energy be- 
havior of the effective impurity problem generated by the 
EDMFT. For this purpose, we choose trial forms for the 
Weiss fields such that the density of states of the fermion 
bath near the Fermi energy is a nonzero constant, 



En 



Nn 



(14) 



while the spectral function of the fluctuating magnetic 
field has a sublincar power-law dependence on energy at 
sufficiently low energies: 



p 



5{uu-Wp) = (Kt/v)\uj\' 



for 



< A, 



(15) 



with 7 < 1. Equivalently — through Eq. (8) — the imagi- 
nary part of the Weiss field takes the form 



where 



Imxo + *0 + ) = C|o>| 7 sgnw, 



C = (K Q g) 2 . 



(16) 



(17) 
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Eq. (15) also defines the cutoff parameter A. The self- 
consistency of our trial forms for the Weiss fields will be 
established in Section IV, where the solution for 7, C, A, 
and No will also be given. 

Having specified the parameters of the impurity model, 
we now proceed to solve this model. Our strategy is first 
to construct the RG equations for the coupling constants 
Jk and g. We then show that a phase transition ex- 
ists as the ratio g/T K is varied. [Here we parametrize 
the Kondo coupling Jk using the single-impurity Kondo 
temperature, 



->-■ 



f-electron 

conduction 
electron 



v/\/\/\/\/\/V vector-boson 

FIG. 2. Bare propagators of the impurity problem defined 
by the Hamiltonian (7), with the spin of the local moment 
represented in terms of (pseudo-)/-electrons. 



1 K 



1 



A)0) 



-1/po(ij.)Jk 



(18) 



where po(fJ-) is the bare conduction-electron density of 
states at the chemical potential The critical value of 
g/T K is determined to linear order in 1 —7. We then cal- 
culate the correlation functions at the critical point, also 
to linear order in 1 — 7. This (1 — 7)-expansion was first 
introduced to the Bose- Fermi Kondo model by Smith and 
Si 39 [following an earlier (1 — 7)-expansion for a spinless 
Bose-Fermi Kondo model 33 ] and by Sengupta. 40 It was 
subsequently used to study a Bose Kondo model with an 
interacting bath by Vojta, Buragohain, and Sachdev. 41 

Both the construction of the RG equations and the 
calculation of correlation functions require a proper han- 
dling of the spin operators. Following Smith and Si, 39 we 
adopt the Abrikosov representation of the spin in terms 
of pseudo-/-electrons, 42 



(19) 



where t x ' v ' z are the Pauli matrices. To stay within the 
Hilbert space for the spin, we assign an energy A to the 
/-electron level and calculate any correlation function 
{V) involving the spin from the corresponding correla- 
tion function in the /-electron representation (D) using 



V = 



Here the prefactor 



lim - el )x V. 

A — *oo 2 
„/3A 



(20) 



is introduced to compensate the 
Boltzmann weighting factor for the singly occupied sub- 
space in the pseudo-fermion basis. Wick's theorem ap- 
plies to T> since it involves only canonical fcrmion opera- 
tors. 

Fig. 2 shows the diagrammatic representation of the 
bare propagators: 



G" f (iu n ) = - 



1 



for the /-electron, 

G b (p, iui n ) 



iuj n — A 



iu> n — E p 



(21) 



(22) 




a ' a ■ 



<7 



a 



FIG. 3. (a) Kondo coupling between the spin of an 
/-electron and the spin of a conduction electron; (b) coupling 
between the spin of an /-electron and the vector-boson field. 



G h Ap,iv„ 



2w v 



(iv m ) 2 - w% 



(23) 



for the vector bosons. G^(p, iv m ) is the Fourier transform 

of G>,r) = -(T T ($ F + $l p )(r)-($ p + $l p )(0))/3. 

In the pseudo-fermion representation, the Hamiltonian 
terms describing the Kondo interaction and the local 
moment-vector boson interaction take the forms 

<7<j' 

(24) 

H 9 = f £^4/^ 2 + ^f\h<t>~ + h.c), 



for the conduction electrons, and 



where a, a' = ±1, $ = £ p (<?p + <£l p ), and 0± = (<j) x ± 
4> V )/V2. Hj and H g are shown graphically in Fig. 3. 



A. RG equations 

We first determine the bare scaling dimensions of the 
couplings Jk and g. S(t) is dimensionless because when 
Jk =5 = 0, (S(t) • S(0))o is a constant. </>(r) has the 
scaling dimension of t~( 1+7 " 2 , and s c (t) that of r _1 . 
Therefore, Jk is marginal, while g has scaling dimension 
(1 — 7)/2, i.e., g is marginal for 7 = 1 and relevant for 
7 < 1. 

We now calculate the (3 function to cubic order in 
the coupling constants. To do so, we first calculate the 
perturbative contributions to the /-electron self-energy, 
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FIG. 4. Diagrams for the /-electron self-energy. 



and the perturbative corrections to the vertices shown in 
Fig. 3. (The self-energies of the fermionic and bosonic 
baths are of order 1/N s n e and vanish in the thermody- 
namic limit for the baths or, equivalently, when the bath 
spectra are continuous.) Details of the calculation are 
given in Appendix A. 

Consider first the /-electron self-energy. The relevant 
diagrams for the perturbative corrections are given in 
Fig. 4. The result is: 



S/(w + A) 



\(N n J K f + \{K,gf 



W 

win—, (25) 

OJ 



where W is a running cutoff energy. Introducing Gf 
through 



Gf(oj + A) = G h f {u + X)Qf((j), 



we have 

g f (u) = i 



- 8 (N a J K f + \{K gf 



(26) 



W 

In—. (27) 

OJ 



Next consider the corrections to the Kondo coupling 
Jk- Due to spin- rotational invariance it suffices to ex- 
amine just the spin-flip processes, say. Also, we set the 
energy of the external /-electrons to w + X and that of the 
external conduction electrons to fi. The diagrams that 
lead to singular corrections to the full vertex Tj(u) are 
given in Fig. 5. We write the full vertex as 



and find that 

7jM = 1 + 



F/M = Jk 7 jM, (28) 



N a J K -]-(N a J K ) 2 -](Kog) 2 In—. (29) 

o 4 OJ 



Finally, consider the corrections to the coupling con- 
stant g. The diagrams contributing singular corrections 
to the corresponding full vertex T g are given in Fig. 6. 
Writing 





FIG. 5. Perturbative corrections to the Kondo coupling 
vertex. 




> 

7f 



FIG. 6. Perturbative corrections to the vertex describing 
the coupling of the local moment to the vector bosons. 



we have 
7sM 



r ff (w) =g 7 fl (w), 



l(N a J K ) 2 + \{K gf 



(30) 



W 

In—. (31) 



We are now in a position to construct the RG equa- 
tions. Lowering the cutoff energy from W to W pro- 
duces the coupling-constant and wave- function renormal- 
izations 

g f (oj, W',J' K ,g') = z f g f (oj, W, J K ,g), 

jj(uj, W\ J' K , g') = Zj^j(oj, W, J Kl g), (32) 

j g (oj,W',J' K ,g') = Zg^gfaW^jcg), 

where 

J' K = zj 1 z J J K , g' = zj 1 z g g. (33) 
This leads to the RG equations 
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FIG. 7. Schematic RG flow of the Bose-Fermi Kondo 
model. Note the unstable fixed point at finite Jk and g. 



J. 



dl 
dg 

-Jl = 9 



K 



N J K 



2 7 2 



-MJ 



2„2 



°K 



_ AT 2 J 2 _ K 2 n 2 



(34) 



where I = ]n(W/W). 

To linear order in 1 — 7, the (3 functions given in 
Eq. (34) yield an unstable fixed point located in the Jk~9 
plane at 



{K g*Y = N J K = {1-^/2. 



(35) 



A separatrix passes through this critical point and 
through the origin, Jk = 9 = 0: to its left (right) in 
the Jk-Q plane the RG flow of Jk is towards infinity 
(zero) (see Fig. 7). For small Jk and g, the position of 
the separatrix, expressed in the form g c (JK), can be de- 
termined by iterating the RG flow. This calculation, the 
details of which appear in Appendix B, gives 



(Xo3c) 2 = 2(l- 7 )exp - 



1-7 
N J< 



K 



(36) 



B. Correlation functions at the critical point 

Consider first the local spin autocorrelation function 
Xloc (r) = (T t S*(t)S*(0)) - lim ± e? x x(r), (37) 



A— -co 2 



where 



x(r) = hT T fl(T)Mr)fl (0)/;(0)), (38) 



which can be calculated using diagrammatic techniques. 
The correlation functions along the entire separatrix 
specified by Eq. (36) are asymptotically the same as their 
counterparts at the fixed point described by Eq. (35). 



As usual, given that 1 — 7 is treated as an expansion 
parameter, the correlation functions at the critical point 
can be calculated via perturbative expansion in the cou- 
plings Jk and g, at the end replacing these couplings by 
their fixed-point values, J* K and g* , respectively. The 
low-order diagrams are shown in Figs. 8(a)-(g). 

Because J* K is of order 1 — 7, the corrections to \\oc 
from the Kondo coupling Jk [Figs. 8(e), 8(f), and 8(g)] 
are of higher-than-linear order in 1 — 7. The calculation 
of the remaining contributions is given in detail in Ap- 
pendix C. Diagram 8(a) contributes to zeroth order in 
I-7: 



X^W = l/4. 



(39) 



Three diagrams contribute to first order in 1 — 7. The 
singular correction from diagram 8(b) is 



(40) 



(Here we consider < r < (5.) The singular contributions 
from diagrams 8(c) and 8(d) sum to 



M)/ x 3 2 sin(7rr//3) 

x[J{r) = --{K g) In —j^. 



(41) 



Collecting all these contributions, noting that (Kog*) 2 = 
(1 — 7)/2, and following the procedure of the standard 
e-expansion, 43 we end up with the following form for the 
local spin correlation function at the critical point: 



Xioc(r) = - 



tt/A/3 



1-7 



(42) 



4 \sin(7rr//3), 
In the zero-temperature limit, this gives, for < 7 < 1, 
1 



xiocM 



r( 7 ) S in^lL_2l(_j W )-7 ) 



2A i_ 7 v,/ - 2 
where T is the gamma function, and for 7 = 

( ^ 1 1 A 
XiocH = — In- 



2A —iu 



(43) 



(44) 



Now we turn to the conduction-electron Green's func- 
tion. The leading nonvanishing contribution to the 
conduction-electron self-energy is given in Fig. 9. Since 
J K is of order 1 — 7, we conclude that 



£ = 0((l- 7 ) 2 ), 



(45) 



i.e., the conduction-electron Green's function is unrenor- 
malized to linear order in 1 — 7. Note that this stands 
in marked contrast to the strong-coupling Kondo phase, 
where E has a pole. 44 

To summarize, for a fixed Jk, the impurity model dis- 
plays a QCP at a critical coupling, g c (JK) specified by 
Eq. (36). The static impurity spin susceptibility diverges 
as 9 — > 9c > as illustrated in Fig. 10(a). At the QCP, 
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FIG. 8. Perturbation series for Xioc 5 defined in Eq. (38), 
showing terms entering at zeroth (a), first (b)-(d), and higher 
(e)-(g) order in 1 — 7. 

/ * 
-^-t "-<--- 



FIG. 9. Leading nonvanishing contribution to the conduc- 
tion-electron self-energy. 

the dynamical susceptibility is given by Eq. (43). Asso- 
ciated with this critical point is the vanishing of a local 
energy scale E* oc , reflecting a critical slowing down. This 
is shown schematically in Fig. 10(b). For g < g c , E* oc is 
finite and serves as an effective Fermi energy scale. 

Subsequent to the work reported here, the RG analysis 
of the Bose-Fermi Kondo model has been extended 45 ' 46 
to higher orders in 1— 7. An unstable fixed point has been 
identified not only in the spin-isotropic model considered 
above, but also in the presence of anisotropy. In both 
cases, the exponent of the local spin correlation function 
at the critical point turns out to be 1 — 7 to all orders in 
1 — 7, in agreement with Eq. (42). 

We note that power-law behavior of local spin- 
correlation functions has also been discussed in impu- 
rity models that arise in the context of quantum spin- 
glasses. 47 ' 40 ' 48 ' 49 There, it is associated with a stable fixed 
point (i.e., a phase) of the impurity model. By contrast, 
what we have discussed above is an unstable fixed point 




FIG. 10. (a) Static local spin susceptibility as a function 
°f 9/Tk- F° r 9 < 9c, the susceptibility is finite and the local 
moment is completely quenched. The divergence of the sus- 
ceptibility at g = g c signals that the dynamics of the local 
moment have become critical, (b) The local energy scale _Ej* oc 
as a function of g/T%. E* oc vanishes at g c , reflecting a critical 
slowing down. 

(a critical point). 

IV. SELF-CONSISTENT SOLUTION 

We must now address two key questions: 

• Can the critical point associated with the local 
problem coincide with the magnetic ordering tran- 
sition? 

• What can cause the fluctuating magnetic field to 
have a sub-ohmic spectrum (7 < 1)? 

The answer to both questions lies in the self- 
consistency requirement. The determining factor turns 
out to be the form of the RKKY density of states de- 
fined in Eq. (5). The significance of pi(e) is that it char- 
acterizes the way the RKKY interaction enters the self- 
consistency condition. Indeed, combining Eq. (5) with 
Eq. (3), we can rewrite the first equation in the self- 
consistency condition (13) as 

^w=hi$fc- (46) 

We consider two types of RKKY density of states p/(e): 
The first type increases from the lower band edge (at 
£ = ^q) with a jump, as shown in Fig. 11(a). The second 
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FIG. 11. (a) An RKKY density of states with a jump at 
the lower edge, e = Iq, characteristic of two-dimensional 
magnetic fluctuations; (b) An RKKY density of states with 
a square-root onset at the lower edge, characteristic of 
three-dimensional magnetic fluctuations. 

type increases from the lower band edge in a square-root 
fashion, i.e., as J~t — Iq . This is illustrated in Fig. 11(b). 

The jump and square-root onset at the lower edge are 
characteristic of magnetic fluctuations in two dimensions 
and three dimensions, respectively. This can be seen, for 
instance, by comparing the form of the RKKY density 
of states for nearest-neighbor real-space coupling on a 
(two-dimensional) square lattice, where 7 q = I (cos q x a + 
cos q y a), and a (three-dimensional) cubic lattice, where 
iq = /(cos q x a + cos q y a + cos q z a). 

In the remainder of this section, we describe the quan- 
tum critical behavior for the two types of RKKY density 
of states. The quantum phase transition takes place as 
the ratio of the RKKY interaction to the Kondo coupling, 



i£Qi 
T° K 



(47) 



is tuned through some critical value S c , which labels the 
QCP. We expect 8 C , the ratio of energy scales charac- 
terizing the competing interactions that give rise to the 
transition, to be of order unity. 

The QCP is signaled by the divergence of the static 
spin susceptibility at the peak wavevector Q. The peak 
susceptibility is related to the spin self-energy M(u) via 
Eq. (3): 



x(Q,w) 



l 



M(uj) + 7 Q ' 
This implies that at the QCP, 

M(w -> 0) = -Iq. 



(48) 



(49) 



A. Self-consistent solution in two dimensions 

With a jump in the RKKY density of states, the self- 
consistency equation Eq. (46) can be written as 



/q+Ao 



de- 



P/(e) 



de- 



MH+e J Iq+Ao M(u) + e 



(50) 



where Ao is defined such that, for e € {Iq, Iq + Ao), P/(e) 
is approximately equal to its value at the lower edge, 
Pj(Iq). For instance, if pi(e) is flat, then A is equal 
to the separation between the upper and lower edge of 
the flat p/(e). For a generic pi(e), Aq is of the order 
of the RKKY interaction as parametrized by, say, |/q|. 
Given Eq. (49), the first term in Eq. (50) is singular; the 
second term, on the other hand, is nonsingular and can be 
neglected. The self-consistency equation then becomes 



XiocM = P/(/q)1ii 



Ar 



M{uj) + Iq ' 



(51) 



Eq. (51), when combined with Eqs. (48) and (49), leads 
to the conclusion that the local susceptibility XiocO^ — > 0) 
becomes singular at the same value of S where the peak 
susceptibility x(Q> w — * 0) becomes singular. As usual, 
the singularity in x(Q, u>) signals the emergence of a crit- 
ical mode associated with the long-wavelength fluctua- 
tions of the order parameter. However, from Section III B 
[see the discussion concerning Fig. 10(a)], we know that 
the divergence of Xioc(w) signifies the emergence of crit- 
ical but spatially local modes; thus, the fluctuations of 
the local moments are also critical. The corresponding 
energy scale E* oc vanishes [Fig. 10(b)]. 

This reasoning leads to one of our key conclusions con- 
cerning the quantum critical behavior of the Kondo lat- 
tice in two dimensions: as we increase the parameter 5, 
spatially local critical modes emerge simultaneously with 
the usual long-wavelength critical modes. 

We now determine the detailed dynamics at the QCP, 
beginning with the conduction-electron self-energy. Since 
the effective impurity model is at its own critical point, 
our result of the previous section [Eq. (45)] implies that 
S vanishes to linear order in 1 — 7. This shows that 
the fermion-bath density of states, assumed to be No in 
Eq. (14), is indeed a nonzero constant and is, in fact, 
equal to the bare conduction-electron density of states 
[see Eq. (2)] at the chemical potential: 



No = Po(aO- 



(52) 



Let us turn to the dynamical spin susceptibility. Com- 
bining Eq. (51) with Eq. (11) yields 



Xo" 1 M + 



1 



xiocM 



= Af(w) 



'Q 



Aq exp 



XiocM 
Pi(lQ) 



(53) 



Inserting Eqs. (16), (36), and (43) into Eq. (53), and de- 
manding a power-law form for the frequency-dependent 
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part of the spin self-energy M(ui), we reach the following 
self-consistent solution for the parameters introduced in 
Eqs. (15)-(17): 



7 = 0+, A = C/k=-T° k 



K ! 



(54) 



where is defined in Eq. (18). The specific form for 
the spectral function of the vector-bosonic Weiss field is 



Imxo 1 {u + i0 + ) = ttA ( In 



A 



while the local susceptibility is 

Xioc(w + «0+) = ^-ln 



sgnw, 



A 



(55) 



(56) 



The spin self-energy is given in terms of xioc through the 
second equality of Eq. (53). The result is 



M(cu) = -I Q 
where the exponent is 



A Hw/A) c 



a = 



2Apj(/ Q ) 



(57) 



(58) 



The parameter A is determined by the Kondo coupling, 
as shown in Eq. (54). The parameter pi(Iq), introduced 
earlier, is the RKKY density of states at the lower edge 
[Fig. 11(a)]. 

We note that the local susceptibility at the critical 
point, given in Eq. (56), is universal. The exponent a 
for the spin self-energy, on the other hand, depends on 
the product Api(Iq). Since the QCP is reached through 
competition between the RKKY and Kondo interactions, 
we expect that this product is close to unity, resulting in 
an a that is not too far from i . The precise value of this 
product, however, is nonunivcrsal: it depends on which 
point of the phase boundary in the two-dimensional (|/q |- 
1%) parameter space is crossed as the system is tuned 
through the quantum phase transition (see Fig. 12). 

The phase boundary can, in principle, be located 
by inserting Eq. (49) into Eq. (53). Recognizing that 
l/Xi oc (w — > 0) = 0, one obtains the condition 



Xo 1 ^ «) = -*: 



(59) 



In this paper we have carried out an asymptotic low- 
energy analysis of our dynamical equations. By contrast, 
Xq 1 (w — > 0) depends (through the Kramers-Kronig rela- 
tion) on the solution at all energies. Therefore, the deter- 
mination of the phase boundary and, hence, the value of 
the exponent a, requires a complete (numerical) treat- 
ment of the problem. Recently, such a treatment has 
been carried out for a Kondo lattice model with Ising 
RKKY couplings. 50,51 A locally critical quantum phase 
transition was identified, 50 with a near-universal expo- 
nent a = 0.72 ± 0.01. Extension of such numerical work 




FIG. 12. Phase diagram of the Kondo lattice in the \Iq \-T% 
parameter space. The solid line is the phase boundary. The 
dashed line shows the tuning trajectory through the quantum 
phase transition for a given system. 



to the isotropic case considered in the present paper re- 
mains an important problem for the future. 

Fig. 12 also illustrates the important point that, in a 
given system, there is only one parameter that must be 
tuned to drive the system through criticality. The dashed 
line represents such a tuning trajectory. 

Two additional remarks are in order. First, the self- 
consistent Weiss field Xo i given in Eq. (55), has a log- 
arithmic dependence. In (imaginary) time, this cor- 
responds to an interaction with an algebraic compo- 
nent, 1/t, multiplied by a logarithmic correction. On 
the other hand, the self-consistent solution for the spin- 
spin correlation function at the critical point, given by 
Eq. (56), corresponds to 1/t, with no logarithmic cor- 
rection. Within the (1 — 7)-expansion carried out in Sec- 
tion III, the logarithmic corrections to the range of the 
interaction do not affect the critical exponent for xioc, 
which remains algebraic without any logarithmic correc- 
tions. We believe that this remains true beyond the 
(1 — 7)-expansion, since the spin-spin correlation func- 
tions at a critical point are in general algebraic and ad- 
mit no logarithmic corrections. For the particular form 
of the self-consistent Weiss-field Xo 1 used here, we can 
demonstrate this point explicitly in a spherical model 
with a long-range interaction specified by x^ 1 ( see Ap- 
pendix D). 

Second, the self-consistency equation (53) dictates 
that, at the critical point, the imaginary part of the 
spin-self-energy ImM(w) is sub-leading compared to both 
Imxo(w) and Imxioc(^)- In principle, we should allow a 
sub- leading term in Xo(^), and determine both the lead- 
ing and the sub- leading terms for the critical xioc(w). 
Calculating the sub-leading terms of xioc(w) is, however, 
beyond the scope of the (1— 7)-expansion. What we have 
done, instead, is to take advantage of the freedom to 
choose the Weiss field such that the self-consistency equa- 
tion is satisfied; Eq. (53), then, allows us to determine 
M from the leading term of xioc(^)- 
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B. Self-consistent solution in three dimensions 

If the RKKY density of states increases at the lower 
band edge in a square-root fashion, 



lim pj(e) oc - Iq, 



(60) 



the self-consistency equation (46) can again be written as 
in Eq. (50), the difference being that A is now defined 
such that the square-root form applies for e € (Iq , Iq + 
A ). Again, for a generic pi(e), A is of the order of the 
RKKY interaction |7q|. 

The second term on the right-hand side of Eq. (50) is 
finite, just as in two dimensions. In three dimensions, 
however, the first term — when combined with Eqs. (49) 
and (60) — also yields a finite value: 



i-Iq 
JIq 



/q+Ao 



(k 



P/(e) 



M(u) + e 



(61) 



Thus, self-consistency dictates that in the three- 
dimensional case, the local susceptibility xioc is finite at 
the QCP. Hence, from Section IIIB [see the discussion 
concerning Fig. 10(a)], we know that the dynamics of 
the local moment are not yet critical. Equivalently, the 
local energy scale E* oc is finite at the QCP. 

We now proceed to extract the critical behavior. To 
be concrete, we work with a semi-circular form for the 
RKKY density of states: 



(62) 



for |e| < I. Solving the self-consistency equation (46) 
yields 



2/ Xioc(w) 



(63) 



which, when combined with Eq. (11), gives 

Xo V) - (//2) 2 Xioc(^). (64) 

The local physics is noncritical at the QCP, so the 
imaginary part of the spin self-energy should reflect the 
decay of spin fluctuations into particle-hole excitations 
and, hence, should be linear in frequency. We can quan- 
tify E* using the definition 28 



ImM _1 (w) = uj/(E* oc ) 2 for \lu\ < E. 



loc ' 



(65) 



Combined with Eq. (49), which now reads M(u^§) 
I, this implies that 



M(w) = / - i{I/El c fu + 0{lo 2 ). 



(66) 



The precise value of E* oc can be determined using a tech- 
nique that is applicable to the strong-coupling regime of 
the Bose- Fermi Kondo Hamiltonian, Eq. (7); one such 
method is slave-boson mean-field theory. 44 



Combining Eqs. (63) and (66) leads to the following 
form for the local susceptibility: 



Xk 



loc 



-2iw/(E? oc y + 0(uj- 3 / 2 ). 



(67) 



Using Eq. (64), the Weiss field is 



Xo V + *0+) 



El 



3/2 



-2i(I/El c fu + 0{^ 2 ). 



(68) 



We should stress that Eqs. (67) and (68) are compati- 
ble. In other words, for a given Bose-Fermi Kondo Hamil- 
tonian [Eq. (7)] and given bath dispersions [Eqs. (14) 
and (15) with 7 = 1/2], the solution will yield an im- 
purity dynamical spin susceptibility whose imaginary 
part is proportional to |w| 1 / 2 sgn(jj. While the correla- 
tion functions on the strong-coupling side of the Bose- 
Fermi Kondo model are not accessible within the (1— 7)- 
expansion (and are in general difficult to determine ex- 
cept for certain large-N limits), we note that the above 
is consistent with the bound set by Griffiths' theorem. 52 
The correlation functions in the paramagnetic phase of 
the one-dimensional long-ranged spherical model provide 
additional support for this conclusion; see Appendix D. 



V. DYNAMICAL SPIN SUSCEPTIBILITY OF 
THE LOCALLY CRITICAL POINT AT FINITE 
TEMPERATURES 

We now return to the locally critical point that occurs 
in the two-dimensional case. Inserting our T = solution 
for the spin self-energy, Eq. (57), into Eq. (3), we arrive 
at the following result for the momentum-dependent dy- 
namical spin susceptibility at T = 0: 



x(q,w,T: 



0) = 



/q+AoHw/A)«' 



(69) 



where A is defined right after Eq. (50). We emphasize 
once again that the local susceptibility, given in Eq. (56), 
is universal. On the other hand, the exponent a is not 
universal, as discussed after Eq. (58). It should also be 
noted that, for q ~ 0, Eq. (69) is expected to be valid 
only if the total spin is not conserved; in heavy fermions, 
this is the case as a result of the strong spin-orbit cou- 
plings. In the conserved case, it would be important to 
keep track of the momentum-dependence of the spin self- 
energy for small q. 

At finite temperatures, the self-consistent solution still 
has 7 = + . From Eq. (42), we have 



Xioc(r) 



I tt/A/3 
4sin(7rr//3)' 



(70) 



yielding 
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Xioc(^ + iO+,T) = 



2A 



In 



A 



2irT 



(71) 



where ip is the digamma function. 

The corresponding spin self-energy can be found using 
the second equality of Eq. (53). The result is 

M (w + iO + ,T) = -7 Q + A a (2TT/A) a T a M(cu/T), (72) 

where 

M(lu/T) = exp[aV(V 2 - iu/2irT)]. (73) 

The momentum-dependent dynamical spin susceptibil- 
ity at finite temperatures follows from inserting Eq. (72) 
into Eq. (3): 



x(q,w,T) = 



Iq + A (2Tr/A) a T a M(Lu/T) ' 



(74) 



lu/T scaling. 26,21 ' 54 ' 55 (This is very different from QCPs 
in insulators in two dimensions, 56,57 where the dynamic 
exponent is usually z = 1 and the corresponding 4> 4 - 
theory is below its upper critical dimension.) 

It should be noted that for two-dimensional magnetic 
fluctuations, the local (i.e., q-averaged) dynamical spin 
susceptibility associated with x G (q, uj) is, in fact, singu- 
lar: 



(77) 



XiocM = dux (q,w) ~ In- 



-IL) 



Within the SDW description [Eq. (75)], the nonlinear 
terms couple only the long-wavelength modes, so this 
singular x? oc IS inconsequential. However, when local de- 
grees of freedom play a central role from the outset — as 
is the case for the local moments in heavy fermions — this 
singularity can cause important nonlinear effects in the 
dynamics of the local objects. The preceding sections 
amount to a microscopic analysis of these nonlinear ef- 
fects. 



VI. GINZBURG-LANDAU THEORY 

A. Gaussian critical point and its instability 

In the picture of the T = SDW transition, there 
is only one type of critical degrees of freedom, namely, 
the long-wavelength fluctuations of the magnetic order 
parameter. The Ginzburg-Landau (GL) action is an ex- 
tension of the standard cf> 4 theory for a classical phase 
transition to allow for order-parameter fluctuations that 
take place not only in space but also in (imaginary) 
time. ' ' 2 ' 53 The temporal fluctuations can be thought 
of as adding z dimensions, where z is the dynamic expo- 
nent. In the antiferromagnetic case (Q being nonzero), 
z = 2 and the action takes the form 

Ssdw= 5iw [ m(q ~ Q, u) ] 

= JdqJ ^[r+(q-Q) 2 + M] [m(q,c)] 2 

+ u n J dcji s ( yi q * 




+... 



(75) 



The effective dimensionality is d c g = d + z. For d = 2 
or d = 3, deft equals or exceeds the upper critical dimen- 
sion of 4, so the critical point is Gaussian. The T = dy- 
namical spin susceptibility assumes the mean-field (RPA) 
form: 



x G (q,^) 



l 



(q-Q) 2 



(76) 



where the linear-in-frequency dependence comes from 
Landau damping (decay into particle-hole pairs). All the 
nonlinear couplings are irrelevant, so the corresponding 
relaxation rate is superlinear, leading to a violation of 



B. Ginzburg-Landau action for the locally critical 
point 

We now address the important question of whether the 
locally critical point is internally consistent beyond the 
microscopic (EDMFT) approximations. 

The defining feature of a locally critical point is the 
coexistence of local and long-wavelength critical modes. 
The long-wavelength modes, as usual, capture the diver- 
gence of the correlation length in space. They reflect the 
fact that there is a broken symmetry in the antiferromag- 
netic metal, and the order parameter corresponding to 
this broken symmetry vanishes continuously as the criti- 
cal point is approached from the ordered side. The local 
critical modes are additional critical modes. They arise 
from a subtle "ordering" on the paramagnetic metal side, 
which also disappears continuously as the critical point is 
approached from the paramagnetic side. The "ordering" 
is associated with the formation of Kondo singlets, a pro- 
cess which yields Kondo resonances and makes the local 
moments part of the electron fluid. It is reflected in the 
formation of a "large" Fermi surface, which encloses a 
volume that counts both the number of conduction elec- 
trons and the number of local moments. 44 The strength 
of these Kondo resonances is captured by a quasiparti- 
cle residue Z, associated with the states on such a large 
Fermi surface. Unlike the SDW case, at the locally criti- 
cal point Z goes to zero over the entire large Fermi sur- 
face. 

The GL action, then, has to be the sum of three parts: 



Slcp = 5 iw [ m (q ~ Q, ] + Sioc + s n 



(78) 



Here, <Si w is the same as in Eq. (75), S\ oc describes the 
couplings among the local critical modes alone, and <S m i x 
takes into account the coupling between the local and 
spatially extended degrees of freedom. 



12 



The notion that local fluctuations are part of the 
quantum critical theory is at first sight surprising, espe- 
cially viewed from what is known about classical (finite- 
temperature) phase transitions. The universality class 
of a classical critical point is entirely determined by 
statics. 58 (Real-time dynamics are considered only af- 
ter the statics have determined the universality class. 59 ) 
Static fluctuations can be energetically favorable only if 
they occur at long distances. The universality class of 
a QCP, by contrast, depends on statics and dynamics 
(quantum fluctuations) on an equal footing. The mixing 
of statics and dynamics opens the door for local fluctu- 
ations to have low energies, provided that their fluctua- 
tions in (imaginary) time are slow. 

We are now in a position to address the robustness 
of the locally critical point. In the microscopic analysis, 
we assumed that the spin self-energy and the electron 
self-energy are independent of the wavevectors q and k. 
Allowing a k dependence in the conduction-electron self- 
energy will not affect our analysis, since we are away 
from half-filling and the conduction-electron density of 
states is expected to remain finite. We now argue that 
allowing a q dependence in the spin self-energy leaves the 
locally critical point intact, provided that a < 1 (as is 
the case observed in CeCu6-xAua; and YbRli2Si2). Con- 
sider first the contribution to the spin self-energy from 
the nonlinear couplings within <Si w . The effective dimen- 
sion d e ff = d + z = 2 + 2/a still exceeds 4, so all the 
nonlinear couplings within S\ w must be irrelevant in the 
RG sense. The q dependence of the spin self-energy will 
then be at most (q— Q) 2 . Consider next the contribution 
to the spin self-energy from S mlK . Since the spin modes 
are coupled to the local modes, the corresponding con- 
tribution cannot have a singular q dependence, either. 
As a result, the spatial anomalous dimension rj takes the 
value 77 = 0, and the zero-temperature, zero-frequency 
spin susceptibility for q ~ Q goes as 



x(q) 



1 



(q-Q) 2 ' 



(79) 



The corresponding local susceptibility remains singular. 
Therefore, the local criticality is robust. 

We emphasize that the locally critical point is non- 
Gaussian because the appropriate field theory is not just 
the usual 4> 4 theory, but instead has the form given in 
Eq. (78) above. Each of the three parts making up Slcp 
contains nonlinear couplings which flow under the RG 
transformation. S\ w (the <p 4 theory) contains nonlinear 
couplings purely between long-wavelength modes. Since 
the effective dimensionality satisfies d a g = d + z > 4, all 
these couplings are irrelevant, and they will contribute 
to the spin relaxation rate a term that is superlinear in 
rp 21,26 However, <Si oc and <S m i x contain nonlinear cou- 
plings involving the local modes. Some of these couplings 
are relevant, and will generate a lincar-in-T contribution 
to the spin relaxation rate. The total relaxation rate 
will be dominated by this second type of contribution, 
and will be linear in T, giving rise to uj/T scaling in 



the dynamical spin susceptibility. Our microscopic the- 
ory provides a specific prescription for determining this 
relaxation rate, as well as the universal scaling functions. 

We close this subsection with a brief comment on the 
general validity of the EDMFT analysis of quantum crit- 
ical behavior in metals. We have seen that, for the 
EDMFT results to be valid, it is crucial 27 that the critical 
point has a mean-field exponent for spatial fluctuations, 
i.e., n = 0. (The temporal fluctuations are allowed to 
have an anomalous exponent.) This requires the non- 
linear couplings among the long-wavelength modes to be 
irrelevant, a condition that is relatively easy to satisfy 
at a T — transition in a metallic system. The results 
presented above serve as a concrete example where the 
spatial fluctuations are Gaussian but the temporal fluctu- 
ations are non-Gaussian. For critical points where r\ 7^ 0, 
by contrast, the EDMFT is not expected to describe 
the correct critical properties. For instance, at finite- 
temperature transitions (where 77 is in general nonzero), 
the EDMFT yields a first-order transition. 60 ' 61 ' 51 ' 62 ' 63 
However, quantum critical dynamics — including those at 
finite temperatures — are controlled by the QCP at T = 0, 
and are correctly captured by EDMFT calculations. 



VII. COMPARISON WITH EXPERIMENTS 

We now compare our conclusions with existing exper- 
iments in hcavy-fermion metals, and also make several 
predictions that can be tested in the future. 



A. Inelastic neutron scattering, magnetization, and 
NMR relaxation rate 

Our theoretical result for the dynamical spin suscepti- 
bility near a locally critical point [Eq. (74)] can be rewrit- 
ten 



x(q,w,T) 



1 



f(d) + AT<*M(lo/TY 



(80) 



where A = Ao(27r/A) Q , and /(q) = I q — Iq vanishes 
at the peak wavevector (q = Q) but is finite at other 
wavevectors. Equation (80) implies that the dynami- 
cal spin susceptibility at the peak wavevector satisfies 
anw/T scaling: 



x(Q,w,T) 



1 



AT a M{u}/ry 



(81) 



where the scaling function A4(uj/T) is given in Eq. (73). 
Eq. (81) is consistent with the inelastic neutron scat- 
tering results 16,17 on CeCue-^Au^ at the critical con- 
centration x c ss 0.1, where it is found that the expo- 
nent a 0.75. (w/T scaling was reported 14 earlier in 
UCus-^Pda; for x = 1 and 1.5. This system, however, 
seems to contain strong disorder, the effects of which re- 
main a subject of debate. 14 ' 15 ' 64 ' 65 ) The scaling function 
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used to fit the experiments is (1 — iuj/aT) a , which be- 
haves very similarly to Eq. (73). The neutron scattering 
results 16 ' 17 at generic wavevectors have also been fitted 
to the form of Eq. (80), with the same exponent a w 0.75. 

In our theory, the locally critical picture arises if the 
magnetic fluctuations are two-dimensional. Within ex- 
perimental resolution, /(q) goes to zero along lines in 
the three-dimensional Brillouin zone, implying that the 
magnetic fluctuations in the quantum critical regime are 
two-dimensional in real space in the frequency and tem- 
perature ranges that have been studied to date. 18 ' 19 ' 16 
It should be noted, however, that the magnetic order- 
ing for x > x c is three-dimensional. To properly address 
the finite-T magnetic ordering transition, it will be nec- 
essary to incorporate the (small) RKKY coupling in the 
third dimension. 66 Whether dimensional crossover even- 
tually takes place at lower temperatures in the quantum 
critical regime is a question that should be addressed in 
future experiments. We hope that our work will pro- 
vide a stimulus for this kind of (challenging) experi- 
ment. It is also desirable that neutron scattering be 
carried out at a QCP where the magnetic fluctuations 
are genuinely three-dimensional. In this regard, we note 
that a violation of lo/T scaling has been reported 9 in 
Cei_ a; La a ;Ru2Si2. 

It also follows from Eq. (80) that the static and uni- 
form spin susceptibility, \ = M/H 7 has a modified Curie- 
Weiss form: 

X(T) = Q+ l BTa > (82) 

where 9 = J q=0 - Iq and B w A A-"(27r) a 
exp[aV'(l/2)]. Eq. (82) fits well 16 ' 17 the magnetization 
data in CeCu6- x Au x . 

There are indications that the modified Curie- Weiss 
susceptibility describes several other heavy-fermion met- 
als near a QCP. 1 One such material is YbRli2Si2, 67 
which exhibits thermodynamic and transport properties 
very similar to those of CeC\i§- x A\\ Xl and has a layered 
structure that suggests the possibility of (quasi-)two- 
dimensional magnetic fluctuations. These characteristics 
make YbRh 2 Si2 a promising candidate for local critical- 
ity. However, inelastic neutron scattering has yet to be 
performed in this system. 

The local susceptibility given in Eq. (71) leads to a 
non-Korringa NMR spin-lattice relaxation rate, 

<~ constant. (83) 

T\ 

(We have assumed that the hyperfine constant is 
not strongly q-dependent.) This prediction might be 
tested using NMR measurements on the Cu sites in 
CeCu6- x Au x or the Si sites in YbRh 2 Si2- 



B. Fermi surface properties, Hall coefficient and 
other transport properties 

On the paramagnetic side of a magnetic ordering QCP, 
the coherence energy scale {E* oc ) is finite. This signals 44 
the formation of heavy quasiparticles, whose Fermi sur- 
face is "large" in the sense that its volume counts both 
the conduction electrons and the local moments. In the 
case of an SDW transition, E* oc remains finite through 
the QCP, and the Fermi surface remains large on the an- 
tifcrromagnctic side (except for the folding of the Fermi 
surface due to a broken translational symmetry). By 
contrast, E* oc vanishes at a locally critical point, and the 
Fermi surface is "small" on the antiferromagnetic side, 68 
i.e., its volume and topology are such that the single- 
electron excitations come entirely from the conduction 
electrons. Thus, the topology of the Fermi surface dif- 
ferentiates the antiferromagnetic metal phase on the or- 
dered side of the locally critical transition from its coun- 
terpart for an SDW transition. Experimental probes of 
the Fermi surface, such as de Haas- van Alp hen measure- 
ments, both inside the antiferromagnetic metal phase and 
close to the QCP, may be used to distinguish between the 
two types of quantum phase transition. 

As a consequence of the large change in the Fermi- 
surface volume across a locally critical point, the Hall 
coefficient at zero temperature is expected to undergo 
a large jump at the QCP. This provides another poten- 
tial test of the locally critical picture. (A more general 
discussion of the behavior of the Hall coefficient can be 
found in Ref. 54; see also Ref. 69.) 

Finally, we consider the longitudinal resistivity. At an 
antiferromagnetic SDW critical point, only certain hot 
spots of the Fermi surface are affected by the critical 
fluctuations. In the clean limit, then, the resistivity is 
expected to retain a Fermi-liquid T 2 form. 70,31 At a lo- 
cally critical point, on the other hand, the entire Fermi 
surface is affected and the resistivity will have a non- 
Fermi-liquid form. 

The change of the Fermi-surface volume is also ex- 
pected to be manifested in the behavior of the residual 
resistivity as a function of the tuning parameter. This 
is best studied by using pressure or a magnetic field as 
the tuning parameter. (By contrast, tuning by chemi- 
cal doping would also induce a significant change in the 
amount of disorder.) We note that the residual resistivity 
does seem to show dramatic changes across the pressure- 
driven QCP in CeCusAu (Ref. 71) and the field-driven 
QCP in YbRh 2 Si 2 (Ref. 72). 

VIII. CONCLUSIONS AND OUTLOOK 

In this paper, we have identified a new class of quan- 
tum critical point in Kondo lattice systems. We have car- 
ried out an EDMFT analysis of the Kondo lattice model 
and found that, when the magnetic fluctuations are two- 
dimensional, a locally critical point arises. Here, critical 
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local modes coexist with the long- wavelength fluctuations 
of the order parameter. We have also argued that the 
locally critical point is robust beyond our microscopic 
calculations. 

Our results explain the salient features of the experi- 
ments in several heavy-fermion systems. In particular, we 
have made comparisons with the existing experiments in 
CeCug-^Au^ and YbRh 2 Si2- We have also made a num- 
ber of predictions concerning the NMR relaxation rate, 
the Fermi-surface volume, and the Hall coefficient. 

The crucial ingredient of local criticality is the pres- 
ence of spatially local critical modes. This contrasts 
with the standard critical theory, in which the only low- 
energy modes are long- wavelength fluctuations of the or- 
der parameter. For heavy fermions, the emergence of 
such critical local fluctuations can ultimately be traced 
to the fact that, as a result of the microscopic Coulomb 
blockade, local moments are "pre-formed" at intermedi- 
ate energies. 73 At an SDW transition, the local moments 
disappear at sufficiently low energies; their quenching by 
the conduction-electron spins leaves only electron-like ex- 
citations (Kondo resonances). At a locally critical point, 
on the other hand, vestiges of the local moments per- 
sist to asymptotically low energies. Such local-moment 
physics is neglected in the SDW treatment of magnetic 
quantum critical points. 

For Mott-Hubbard systems, the microscopic Coulomb 
blockade is responsible for the formation of the Mott 
insulator. One is therefore led to speculate that low- 
energy local modes may well play a key role in met- 
als near a Mott transition. This issue is becoming in- 
creasingly important due to the observations of (nearly) 
quantum critical behavior in transition-metal oxides, 74 as 
well as the apparent scaling behavior in high-temperature 
superconductors . 3 ' 4 ' 75 ' 76 
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Consider first the pseudo- /-electron self-energy. There 
are two contributions, as specified by Fig. 4. The contri- 
bution (a) is given by 

3 j: 



p,p' lu>\ lu>2 

xG b (p, iu) 2 )G b (p', iuif + iuj 2 - iui)- (Al) 

Here the factor of 3 comes from summation over the three 
spin components. Carrying through the summation over 
iu>i and zw2, setting iuif = lo + A, and using the fact that 
A — > oo, we have 

^\u+\)=lj%j2 f{£p)[1 ~ f{ei " )] 



lj + e p - e p > 



= ~(N J K f 



D 



2D In 2 + ui In — + OUJ) 

LO 



, (A2) 



where D is the half-bandwidth, and the Fermi factor 
/(e) = (e^ 6 + Note that in the second line we 

have set T = 0, since our purpose here is to construct 
the RG equations. 
Similarly, 

^N^nEE^N^^-^). (A3) 

p iui 

To linear order in 1 — 7, we can set 7 = 1 in and the 
corresponding spectral function is 



Me) 



--ImVG^,e + iO+) 

P 

K$e for |e| < A. 



(A4) 



We then have 



£< 6 >(w + A)=| 5 2 / deA^e)- 



n B (e) 



A + wln— + 0(lu) 



(A5) 



where the Bose factor ns(e) = (e /3e — The constant 

parts can be absorbed in a shift of A. The lo lnw parts are 
retained in Eq. (25), where we have replaced the cutoffs 
by the running cutoff W . 

We now turn to the corrections to the Kondo coupling, 
as given in Fig. 5. We set the external frequencies of the 
conduction electrons to 0, and those for the pseudo- /- 
electrons to iuif. The contribution from Fig. 5(a) is 



APPENDIX A: DERIVATION OF THE RG 
EQUATIONS 

In this appendix, we present the derivation of the RG 
equations (34), following a procedure outlined in Ref. 39. 



LO + € 2 UJ 



(A6) 



where, in the second equality (and at similar places in the 
equations to follow) we take A — > 00. The contribution 
from Fig. 5(b) is 
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r>( fe )c i ^ Jk s^bf \nbf ■ \ i From Eq. (33), we have 



2/3 f 



& - - — = -J^JVbln-. (A7) 



The contribution from Fig. 5(c) is 



Jjr — Jk ~\~ Jk 



1 



n J k - MJ K y - {K o9 y 



In- 



.9 = .9 -.9 



-(7VoJx) 2 + (^og)' 



In 



W_ 

w 



w 
W 

(A16) 



T^ c) ^ f )= -^J2xc(^f -^[G^)] 2 , (A8) 

where Xc is the bare local conduction-electron sus- 
ceptibility, which has a spectral function 



The result is 



= l x »(e + iO+)=N 2 e. 

7T 



(A9) 



l+n B {e) 



(A10) 



rW( W + A)=-i4 JdeA Xc (e) 

1 , , D 
= -8^ In-. 

The contribution from Fig. 5(d) is 

r ? ) N)=irEE G J(P' iw /- iw i)[ G /^)] 2 - ( An ) 

Using Eq. (A4), we have 



T { f\u + \)=-\g 2 J K J deA^e) 



l+n B (e) 



-]j K (K g) 2 In-. 

4 tt> 



(A12) 



Finally, we consider the corrections to the coupling 
constant g. We set the external frequency for the <j) prop- 
agator to 0, and that for each pseudo-/-electron propa- 
gator to iojf. Then 

r W(o;+A)=-^J^Xc(^ / -ia; 1 )^(ia; 1 )] 2 | ia , /=w+A 



= -^?(AW 2 ln-, 



(A13) 



and 



r( b )(o, + A)=|^^G^^ / - J c 1 )[G^^ 1 )] 2 |^_ +A 

(A14) 



\g{K g) 2 \n-. 
4 w 



Setting oj = W in Eq. (32) leads to 



«/ = 1 + 
zj = 1 + 



Zg = 1 



In 



W_ 
W> 



N J K - l(N J K f - \(K g) 2 



W 

^- (A15) 



I(NoJk) 2 + \{K gf 



i W 



APPENDIX B: DETERMINATION OF THE 
SEPARATRIX IN THE J K -g PLANE FROM THE 
RG FLOW 

Our starting point is the RG equations (34), which 
yield the RG flow shown in Fig. 7. We imagine approach- 
ing the scparatrix from the local-moment side. The sep- 
aratrix is identified when the initial parameters are such 
that NoJk no longer flows to zer o but inst ead is renor- 
malized to its fixed-point value \J {I — 7)/2. 

Consider the separatrix near the origin. Provided that 
we work either on the local-moment side or on the sepa- 
ratrix itself, (N Jk) 2 is at most of the order of (1 — 7) 2 
during the entire flow and can safely be ignored inside 
the brackets on the right-hand sides of both flow equa- 
tions (34). (This is not the case on the Kondo side.) 
With this approximation, and the shorthand notation 



a = N J K , b=(K g) 2 , 
the RG equations simplify to 
da ( u\ 



1 



' 7 



ir 2b 



These equations have an unstable fixed point at 
a* = b* = 1. 



Solving Eq. (B3) yields 
, 1-7 



1 + (c - l)e 



-0-*o)(l-7) 



Here. 



co 



1-7 
2& ' 



(Bl) 

(B2) 
(B3) 

(B4) 

(B5) 

(B6) 



and 60 and Iq (as well as a , z , and y used below) are 
the initial values of the corresponding parameters. 
We can now proceed to solve Eq. (B2). Let us write 



a = zy, 

where the factor z is defined as satisfying 
dz 



dl 



= -bz. 



(B7) 



(B8) 
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When combined with Eq. (B5), Eq. (B8) can be easily 
solved to give 



z = z 



c e 



-(;-; )(i- 7 ) 



(B9) 



l + (c -l)e-( i - ( o)(i-7) 
The other factor, y, satisfies the differential equation 

(B10 » 

which, when combined with Eq. (B9), leads to 

-l 



y = yo 



i 



1-7 



where 



,(I-I )(1_ 7 ) 



^4(0 



dX 



AV A + c - 1 



(Bll) 



(B12) 



Using 1 



,S2 



/ 



dx 



y/x(a + bx) iy/aS 



In 



a — bx + 2iV abx 
a + bx 
for ab < 0, 



(B13) 



and recognizing that, for parameters near the origin, Co = 
(1 — 7)/(26o) <C 1, we obtain the following asymptotic 
form in the limit I — Iq 3> 1/(1 — 7): 

AQ -» 00) w -L ln(4co) - 2 e -^-^ 1 -^/ 2 . (B14) 

Combining Eqs. (B7), (B9), (Bll), and (B14), we arrive 
at the asymptotic solution for a: 



a(l — ► 00) 



QWCoe 



-(l-lo)(l— Y)/2 



5+ a "vf° e -('-'o)(l-7)/2' 



where 



B = 1 



o 



1 



7 



■ln(4co). 



(B15) 



(B16) 



It is now clear that when B is positive, a(l = 00) = 0. In 
other words, the system falls in the local-moment regime. 
On the other hand, a(l = oo) = a* when B = 0, i.e., when 



a ln(4c ) = 1-7. 



(B17) 



This condition specifies the separatrix. Written in terms 
of the original variables, through Eqs. (Bl) and (B6), this 
becomes Eq. (36). 

For completeness, wc also write down the equation for 
the separatrix when the condition c S> 1 is not satisfied. 
In this case, a(l) still has the asymptotic form specified 
by Eq. (B15), but Eq. (B16) is replaced by 



APPENDIX C: SPIN-SPIN CORRELATION 
FUNCTION AT THE CRITICAL POINT: 

(l-7)-EXPANSION 

Here we calculation the spin-spin correlation function 
Xloc (r) ee I<T T S-(r)S+(0)) = lim V*, (CI) 

A A^oo Z 

where \ lSi calculated in the pseudo-fermion representa- 
tion: 



1 



X(r) = ? (T T /|/ T (r).f T t / i (0)). 



(C2) 



It turns out to be convenient to calculate x m terms of 

G b (r) = e"^[(/(A) - l)6(r) + /(A)6(-r)]. (C3) 

The contributions to \ are given in Fig. 8. Note that 
any diagram containing more than one isolated /-electron 
bubble vanishes once the limit A — > 00 is taken. Fig. 8(a) 
yields 



X {a) {r) = - l -G){r)G){-r). 
In the following, we will restrict < r < (i. Then 



X (a) {r) 



(C4) 



(C5) 



The contribution from Fig. 8(b) is 

X {b) (r) = £ [ P dn f dT 2 Gl(r 1+ T 2 ) x 

x G b f ( Tl )G b f (T 2 )G b f (-n-T)G b f (-t 2 +t), (C6) 
and those from Figs. 8(c) and 8(d) sum to 

X (cd) {r) = -^|! jf" dn j* dr 2 G\{r - n - r 2 ) x 

xGjfnjGj^JG^r - n - r 2 )G)(-r) (C7) 

Anticipating that g 2 is of order 1 — 7, to linear order in 
1 — 7 wc can use Eq. (A4) , which corresponds to 

OU^fdeMe^-l-e— 



= Kl 



Vsin(7rr//3) 



over the range 



Irl >> \ and 1/3 - rl >> 4- 
A 



B 



A' 



(C8) 



(C9) 



«oV c o 



In 



2co-l + 2v/co(co-l) , (B18) We then have 



(l- 7 )V^o _ T 

and Eq. (B17) is modified accordingly. It is straightfor- 
ward to check that the choice (ao, bo) = (a*,b*) satisfies 
the condition B = 0. In other words, the unstable fixed 
point is part of the separatrix, as it should be. 



X (fc) W- -^ e ~ PX [ T dn£ dT 2 G b 4T 1+ T 2 ) 



(C10) 
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where the factor A appears due to the fact that Eq. (C8) 
is valid only in the range specified by Eq. (C9). Similarly, 

X M) (r)= ^e"^ jf dn £ T1 dr 2 G%{r - n - r 2 ) 



(Cll) 



APPENDIX D: LONG-RANGED 
ONE-DIMENSIONAL SPHERICAL MODEL 
WITH A LOGARITHMIC-CORRECTION TO ITS 
RANGE 

Consider a spherical model, defined on a chain of M 
sites, with a reduced Hamiltonian 



1~Lsp ~ „ ^""^ VjjSjSj. 



(Dl) 



Here, Si is a scalar variable constrained by the condition 

M 



]>> 2 = M, 



(D2) 



i=l 



with M being taken to infinity at the end of the calcu- 
lation. The coupling Vij represents a long-ranged inter- 
action whose Fourier transform v(p) = Vije~ lpRij is, 
for small p, either 



v(p) = v — vi\p\ x with < x < 1, 



or 



v(p) = v - 



Hi/\ P \y 



(D3) 



(D4) 



in either case, v and v\ are positive constants. 

The long-range interaction corresponding to Eq. (D3) 
has the form 



consistent with those for the paramagnetic side of any 
generic one-dimensional classical spin models with power- 
law long-ranged interactions. 78-80 

We now consider the problem with a long-range inter- 
action specified by Eq. (D4). In real space, v(r) has a 
l/\r\ dependence with a logarithmic correction. We ex- 
pect that, at the critical point, the spin-spin correlation 
function should have a pure power-law decay without any 
logarithmic correction. In the following, we show that 
this is indeed the case. Standard manipulation 77,81 leads 
to the following form for the free-energy: 

m = Xo-ljd P ln(xo-v + ] ^), (D8) 

where the Lagrange multiplier, Ao, introduced to enforce 
the constraint, Eq. (D2), is determined by minimizing F: 



u 

2 J X -v 



dp 



+ VMi/H) 



= i. 



(D9) 



The critical point occurs when the solution is such that 
A = v, corresponding to 



MVbl) _ L 



(D10) 



Note that the integral on the left hand side is infrared 
convergent, establishing the existence of a phase tran- 
sition. The correlation function at the critical point is 
given by 



v{ \p\ 

In real space, this corresponds to 

x cri M ~ A 



(Dll) 



(D12) 



without any logarithmic correction, which is what we set 
out to establish. 



v(r) 



\l+x ' 



(D5) 



The solution to this problem is known. 77 An unstable 
fixed point separates a paramagnetic state from an or- 
dered state. The spin-spin correlation function at the 
critical point behaves as 



1 



\l-x ' 



while that on the paramagnetic side varies as 

1 



x par » 



\l+x ' 



(D6) 



(D7) 



Eq. (D7) implies that on the paramagnetic side, the ex- 
ponent of the spin-spin correlation function saturates 
the bound set by Griffiths' theorem. 52 The result is also 
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